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The scaling theory of Anderson localization is based on a global conductance ql that remains a 
random variable of order O(l) at criticality. One realization of such a conductance is the Landauer 
transmission for many transverse channels. On the other hand, the statistics of the one-channel Lan- 
dauer transmission between two local probes is described by a multifractal spectrum that can be 
related to the singularity spectrum of individual eigenstates. To better understand the relations be- 
tween these two types of results, we consider various scattering geometries that interpolate between 
these two cases and analyse the statistics of the corresponding transmissions. We present detailed 
numerical results for the power-law random banded matrices (PRBM model). Our conclusions are : 
(i) in the presence of one isolated incoming wire and many outgoing wires, the transmission has the 
same multifractal statistics as the local density of states of the site where the incoming wire arrives; 
' (ii) in the presence of backward scattering channels with respect to the case (i), the statistics of the 

, transmission is not multifractal anymore, but becomes monofractal. Finally, we also describe how 

. these scattering geometries influence the statistics of the transmission off criticality. 
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INTRODUCTION 



Within the field of Anderson localization [l[ (see the reviews [U, H, EL S B S HI ) > one main achievement has been 
the formulation of a scaling theory (9J describing the renormalization flow of the typical global conductance as 
a function of the linear size L of the system. The critical point then corresponds to a finite conductance, whereas 
the delocalized phase corresponds to the growing RG flow gi, ~ L d ~ 2 in dimension d > 2, and the localized phase 
corresponds to the decaying RG flow gi, ~ e~^ cst ^ L . Although this scaling theory was first formulated in terms of 
the Thouless definition of conductance, based on the sensitivity to boundary conditions, it was then realized that a 
more straightforward definition of a global conductance can be obtained via the quantum scattering theory where the 
quantum transmission T|™ c ' is g iven by the many-channel (m.c.) Landauer formula [10. IllL Il2l. I IL 31] that generalizes 
O \ the one-dimensional case 0, [HI EH • 

This scattering definition has the advantage to have a well defined meaning for each disordered sample, and thus 
I \ one may study its probability distribution over the samples (lp| (see Q for a review of the results concerning this 
■ distribution both in the localized phase and at criticality) . From the point of view of the scaling theory, the important 
' point is that at criticality, the many-channel Landauer transmission T™" c ' remains a random variable of order O(l) as 
L — > +oo. On the other hand, it is now well established that individual critical eigenfunctions are multifractal (see the 

.' reviews [6, 8]). As a consequence, the one-channel Landauer transmission T^ 1 ' 1 ' between two local probes is expected 
to display also multifractal properties at criticality [l7], EH Ell- To better understand the relation between this 

mnltifrar>tnl ctntictir>a nf tlio nnp-rKflnnAl T .and anor frimcmiCQinn ^(I'-O 



multifractal statistics of the one-channel Landauer transmission T L > and the many-channel Landauer transmission 
T™ c -, we consider in this paper various scattering geometries that interpolate between the two in order to characterize 
the statistics of the transmission in these intermediate conditions. We present detailed numerical results for the 
Anderson critical points of the Power-law Random Banded Matrix (PRBM) model, where the criticality condition is 
. £^ \ exactly known. 

The paper is organized as follows. In section [II] we recall the definition and some properties of the PRBM model. 
We then consider the statistics of the transmission at criticality for the following scattering geometries : 

- one incoming wire and one outgoing wire (section IIIip 

- one incoming wire and many outgoing wires (section IIV() 

- one incoming wire with many backward channels and many outgoing wires (section |V| 

For completeness in section IVI1 we describe how these scattering geometries influence the statistics of the trans- 
mission off criticality. Our conclusions are summarized in section IVIII A brief reminder on multifractality of critical 
eigenfunctions is given in Appendix A. 



II. REMINDER ON THE POWER-LAW RANDOM BANDED MATRIX (PRBM) MODEL 

Beside the usual short-range Anderson tight-binding model in finite dimension d, other models displaying Anderson 
localization have been studied, in particular the Power-law Random Banded Matrix (PRBM) model, which can be 
viewed as a one-dimensional model with long-ranged random hopping decaying as a power-law (b/r) a of the distance 
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r with exponent a and parameter b. Many properties of the Anderson transition at a = 1 between localized states 
with integrable power-law tails (a > 1) and extended states (a < 1) have been studied, in particular the statistics 
of eigenvalues J2ll HI, H3], and the multifractality of eigenfunctions [HI, HI, US H3, 0, Ha], including boundary 
multifractality [301 ]. The critical properties at a = 1 depend continuously on the parameter b, which plays a role analog 
to the dimension d in short-range Anderson transitions [24| : the limit 6^1 corresponds to weak multifractality ( 
analogous to the case d = 2 + e) and can be studied via the mapping onto a non-linear sigma-model [2(| , whereas the 
case & < 1 corresponds to strong multifractality ( analogous to the case of high dimension d) and can be studied via 
Levitov renormalization [24], |32| • Other values of b have been studied numerically [2J, |25|, |26|, |27[ . The statistics of 
scattering phases, Wigner delay times and resonance widths in the presence of one external wire have been discussed 
in [3l|, UK . Related studies concern dynamical aspects [34j , the case with no on-site energies [IH] , and the case of 
power-law hopping terms in dimension d > 1 [H, H3, [H[ . 

In this paper, we consider the PRBM in dimension d = 1 cither in the line geometry with two boundaries or in the 
ring geometry with periodic boundary conditions, in the presence of various external wires to measure the transmission 
properties. 



A. PRBM with boundaries (line geometry) 

In the line geometry, the distance r%j between two sites is simply 

n,j = \i-j\ (1) 

The ensemble of power-law random banded matrices of size L x L is then defined as follows : the matrix elements 
Hi t j are independent Gaussian variables of zero-mean Hij = and of variance 
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B. PRBM with periodic boundary conditions (ring geometry) 



In the ring geometry with periodic boundary conditions, the appropriate distance r^j between the sites i and j on 
a ring of L sites is defined as [24[ 

ry = - S1 n(— _j (3) 

The matrix elements ffjj are then defined as before in terms of this distance (Eq. [5]). The main property of this ring 
geometry is that all sites are equivalent, whereas in the line geometry there are two boundaries at i = 1 and at i = L. 



C. Scattering transmission in the presence of external wires 



In the following, we will consider various scattering geometries, where the disordered sample is linked to one 
incoming wire, parametrized by the half- line (x < x m ), and to one or many outgoing wires, parametrized by the 
half-lines (xj > x° ut ). In each case, we are interested in the eigenstate \ip > that satisfies the Schrodinger equation 

H\^P>=E\^> (4) 

inside the disorder sample characterized by the random -Hjj, and in the perfect wires characterized by no on-site 
energy and by hopping unity between nearest neighbors. Within these perfect wires, one requires the plane-wave 
forms 

i>in(x < x m ) = e ^-* in ) + re -ik(x-x in ) 
Tpout(xj > <) =t 3 e^- x T t ) (5) 

These boundary conditions define the reflection amplitude r of the incoming wire and the transmission amplitudes 
(tj) of the outgoing wires. 
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The Landauer transmission T of a subset J of the outgoing wires defined as 

j&J 

is then a number in the interval [0, 1]. The various scattering geometries considered below differ in the number of the 
possible out-going wires and in the subset J used to measure the transmission. 

To satisfy the Schrodinger Eq. [4] within the wires described by Eq. [5j one has the following relation between the 
energy E and the wave vector k 

E = 2cosk (7) 

To simplify the discussion, we will focus in this paper on the case of zero-energy E = (wave- vector k = tt/2) that 
corresponds to the center of the band. 

In the following, we study numerically the statistical properties of the Landauer transmission T for systems of sizes 
50 < L < 1800 with corresponding statistics of 10.10 8 > n s (L) > 2400 independent samples. For typical values, the 
number n s (L) of samples is sufficient even for the bigger sizes, whereas for the measure of multifractal spectrum, we 
have used only the smaller sizes where the statistics of samples was sufficient to measure correctly the rare events. To 
measure numerically multifractal spectra, we have used the standard method based on g-measures of Ref. [39| (see 
Appendix B of 19] for more details). 



III. STATISTICS OF THE ONE-CHANNEL TRANSMISSION t£ a) 
A. Scattering geometry 

The one-channel transmission T^' 1 ' is the Landauer transmission when the disordered sample is connected to one 

incoming wire and one outgoing wire, see Fig. Q]for the case of the PRBM model. This one-channel transmission T^ 1 ' 1 ^ 
is a very interesting observable to characterize Anderson transitions [UdH : (i) it remains finite in the thermodynamic 
limit only in the delocalized phase, so that it represents an appropriate order parameter for the conducting/non- 
conducting transition; (ii) exactly at criticality, it displays multifractal properties in direct correspondence with the 
multifractality of critical eigenstates, i.e. it displays strong fluctuations that are not captured by more global definitions 
of conductance. 



B. Multifractal statistics of Tj^ at criticality 

As first discussed in pjj for the special case of the two dimensional quantum Hall transition, the critical probability 
distribution of the one-channel transmission T^ 1 is described by a multifractal spectrum <&( 1 » 1 )(/c) 

Prob f t£> 1) ~ L- K ) dT cx L* C1,1) (*>d« (8) 
Its moments involve non-trivial exponents X^ ljl \q) 



(if 1 ')' ~ / dKL* 11 ' 1 ^* ^L-x™^ (9) 
that can be computed via the saddle-point in n 

- A (M) (q) = max [$ (M) (k) - qn] (10) 
K L J 

As stressed in [17j], the physical bound J^ 1 < 1 on the transmission implies that the multifractal spectrum exists 
only for k > 0, and this termination at k = leads to a complete freezing of the moments exponents 

X^(q)=X^(q sat ) for q > q sat (11) 

at the value q sat where the saddle-point of the integral of Eq. [5] vanishes n{q > q sat ) = 0. 
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(a) Line geometry 
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(b) Ring geometry 



,„- e ikx +re ikx 
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ifiout = te 



ikx 



disordered sample 

FIG. 1: Scattering geometries defining the one-channel transmission t£ (a) Line geometry : the incoming wire is attached to 
site 1, and the outgoing wire to site L. These two sites are 'boundary sites' for the disordered sample (they have less neighbors 
than bulk sites), (b) Ring geometry : the incoming wire is attached to site L, and the outgoing wire to site L/2. In the ring 
geometry, all sites are 'bulk sites' for the disordered sample (all sites have the same number of neighbors). 



It is very natural to expect that <I>( 1,1 ^(k) is related to the multifractal spectrum f(a) concerning eigenfunctions 
(see the Appendix [5] for a brief reminder) . The possibility proposed in [17| is that before the freezing of Eq. [11] 
occurs, the transmission should scale as the product of two independent local densities of states (Eq. IA6|) 



X^\q) = 2A(q) for q < q s 



(12) 



We refer to [13] for physical arguments in favor of this relation. Equations [TT] and [T2] for the moments exponents are 
equivalent to following relation between the two multifractal spectra 



^'^(k > 0) = 2 \yf{a = d+~)-d 



In particular, the typical transmission 



rptyp 
1 L 



is expected to decay at criticality with some power-law 



T. 



typ 



x 



JnT L 



(i.i) 



where the exponent niyp^ ^ s the point where <&( 1,:l )(k;) reaches its maximum $( 1 4)( / 
directly related via 



typ j 



(Id) 
Hyp 



(13) 
(14) 
(15) 

From Eq. [H it is 
(16) 
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to the typical exponent at vp that characterizes the typical weight of eigenfunctions 

mtiypOCj^; (17) 

In our recent work , we have tested in detail these predictions for the statistics of the one-channel transmission 
T^ 1 ' 1 ^ for the Power-law Random Banded Matrix (PRBM) model where d = 1 in the ring geometry (see the scattering 
geometry (b) of Figure [T]), where all sites are 'bulk sites'. We thus refer the reader to [19| for detailed numerical data. 




K q 



FIG. 2: Multifractal statistics of the one-channel transmission T^' 1 ^ for b — 0.1 at criticality a = 1 for the line geometry 
and for the ring geometry (see Fig. [TJ (a) the multifractal spectra $^ 1,1 '(/t) corresponding to the Ring geometry and ^surf( K ) 
corresponding to the Line geometry are maximal at the typical values K tvp ~ 1.33 and k*^./ — 1-64. (b) the corresponding 
moments exponents X' 1 ' 1 ' (g) and X^J(q) saturate at the values X^'^fa > q sa t) = — 5^ 1 ' 1 '(0) ~ 0.58 and X^j(q > q sa t) = 
-*<$(0) ~ 0.69 

However, in other localization models it is usual to attach leads to the boundaries of the disordered samples. It is 
thus interesting to consider also the scattering geometry (a) of Figure Q] where the incoming wire and outgoing wire 
are attached at boundary points, since one expects that at criticality, boundary points are characterized by a different 
multifractal spectrum f sur f (a) [4(| IS, EH ■ One then expects that the relation of Eq. [T3] becomes 

> 0) = 2 [f surf (a = d+ |) - d s ] (18) 

when the transition is measured between two boundary points, where d is the bulk dimension and d s is the surface 
dimension. Here the line geometry (a) of Figure [1] corresponds to d = 1 and d s = 0. On Fig. ® we compare for 
the case b = 0.1 at criticality a = 1 the bulk multifractal spectrum ^ 1,1 '(k) and the surface multifractal spectrum 
^sir/W - W e find in particular that the typical exponent for the surface case n[yp S>sur ^ (b — 0.1) ~ 1.64 coincides 

numerically with the typical value of the bulk case {b = 0.05) ~ 1.64 measured in This is in agreement with 
Eq. 29 of Ref. [4l| stating that in the regime of small b, the surface multifractal spectrum of parameter b coincides 
with the bulk multifractal spectrum of parameter 6/2 (here we consider the case p = in the notations of Ref. [4ll|). 
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FIG. 3: Scattering geometry denning the 'one-to-many' transmission T£ for the Ring geometry : the incoming wire is 
attached to site L, and outgoing wires are attached at all sites j satisfying L/4 <j< 3L/4. The total transmission is given by 
Eq. [19] hi terms of the amplitudes tj of the outgoing wires and of the reflexion amplitude r of the incoming wire. 



IV. STATISTICS OF THE ONE-TO-MANY-CHANNEL TRANSMISSION T^ 1 ' 



A. Scattering geometry 



We now consider the scattering geometry of Fig. [3] where the incoming wire is attached to site i = L, and where 
there are (L/2+1) outgoing wires attached to the sites j = L/4, L/A+l, 3L/4. We are interested in the transmission 

3L/4 

T { L hm) = £ \t 3 \ 2 = l-\r\ 2 (19) 
where the tj are the transmission amplitudes and where r is the reflexion amplitude of the incoming wire (Eqs[5]). 



B. Multifractal statistics of T^ 1,m ' at criticality 



In this scattering geometry, one expects that that the critical statistics is still multifractal with a spectrum $( 1 ^ m ) (k) 

-(1.™) _ T.-A AT nr r.* (1 ' m) («) 

and non-trivial exponents X t ^ 1 ' m \q) 



Prob(T™ ~L- K ) dT l oc L [K} dn (20) 



(r^- m) ) ~ j dKL* ll ' m) M-** ^i-x^m (21) 

Again the physical bound T9~' m ' < 1 on the transmission implies that the multifractal spectrum exists only for k > 0, 
and this termination at k = leads to a complete freezing of the moments exponents 

X^ m \q) = X^ m \q sat ) for q > q sat (22) 
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FIG. 4: Multifractal statistics of the one-to-many transmission for b = 0.1 at criticality a = 1 for the scattering geometry of Fig. 
0(a) the multifractal spectrum <J> ( - 1 ' m ) (k) as compared to the multifractal spectrum $ < - 1 ' 1 ' (k). The typical values /4»p™' ~ 0.66 

and n^f ~ 1.32 satisfy the relation of Eq. [26] The termination values $< 1 ' m '(0) ~ -0.27 and $ (M) (0) ~ -0.58 satisfy the 
relation of Eq. 1271 The full relation between the two multifractal spectra is given by Eq. (b) the corresponding moments 
exponents X^ 1,m \q) as compared to X' 1 ' 1 ' (g). 



at the value q sa t where the saddle-point of the integral of Eq. 1211 vanishes K,(q > q sa t) — 0. However in the presence 
of many wires (more precisely, whenever the number of outgoing wires is greater than L* ( ( K=0 )/ 2 according to the 
discussion of the previous section), one expects that the relation of Eq. [T2l becomes 

X^ m \q) = A(q) for q < q sat (23) 

meaning that the transmission simply scales as the local density of states seen by the incoming wire. Equations [H] 
and [221 for the moments exponents are equivalent to following relation with the eigenfunction singularity spectrum 
/(«) 

$( 1,m ) (k > 0) = f(a = d + k) - d (24) 
Equivalently using Eq. 1131 one obtain the simple relation 

$ (^) (K >0) = ^M (25) 
In particular, the typical values of k (where the $ reach their maximal value 0) are related by 

4l P m) = 4^ (*) 

The termination values at k = arc also related by the simple relation 

(0) = (27) 

We have checked these relations (Eqs 12^ 1 126 ) I27 p for the case b = 0.1 at criticality a= 1 (see Fig. [3J. 
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FIG. 5: Scattering geometry defining the 'one-to-many' transmission j^ is ' m > m the presence of backward channels : in addition 
to Fig. [3] there are now outgoing wires for 1 < j < L/4 and 3L/4 < j < L — 1 that are considered as backward scattering 
channels. The total transmission is given by Eq. [28] in terms of the amplitudes tj of the outgoing wires for L/4 < j < 3L/4. 



V. STATISTICS OF THE ONE-TO-MANY-CHANNEL TRANSMISSION T ( L 1B ' m) IN THE PRESENCE 

OF BACKWARD CHANNELS 



A. Scattering geometry 

We now consider the scattering geometry of Fig. [5] : the only difference with the previous case of Fig. [3] is the 
presence of backward scattering channels near the incoming wire. We are interested into the transmission 

3L/4 

T £B, m) = J2 \ tj \* (28) 

3=L/4 

where the tj are the transmission amplitudes (Eqs [5]) . 



B. Non-multifractal statistics at criticality 



In this case, we find numerically that the statistics of the transmission is not multifractal anymore, but mono-fractal 

T (lB,m) _ (2Q) 

where r is a random variable of order 0(1). Moreover, the monofractal exponent is 

Kback(b) = 1 for all b (30) 

as shown on Fig. [S] for the three values b — 0.01, b = 0.1 and b = 0.25, i.e. it does not depend anymore on 
the multifractal critical properties that are known to change continuously with the parameter b in the PRBM. For 
instance for the one-channel transmission, we have measured in [l9| the typical exponents K&p(b = 0.01) = 1.92, 
n[yp\b = 0.1) = 1.33, and K^j/p 0> = 0.25) = 0.77 which are significantly different. 



In L 



FIG. 6: Statistics of the transmission j£ 1S for the scattering geometry of Fig. \5\ at criticality a = 1 : lnT^ = lnTz, as a 
function of InL. The slopes for the three values b = 0.01, b = 0.1 and b — 0.25 coincide Kbackib) = 1 (see Eq. [30}. 



C. Consequences for the usual many-channel (m.c.) transmission T™' c ' 

For the short-range Anderson tight-binding model in dimension d — 3, the usual many-channel Landauer transmis- 
sion T™ c ' corresponds to the incoherent sum 

T? °- = £ T * (31) 

i£L d s 

of L rfs = contributions Ti, where Ti is the coherent scattering transmission for one- incoming wire arriving at 

i, with L ds backward scattering channels and with L ds forward scattering channels (see for instance Fig. 2 of Ref 
[l3|). The statistics of this many-channel transmission has been much studied (see the review 0) : T™ c - is a random 
variable of order 0(1). One of our motivations for the present study was to determine whether the statistics of each 
contribution Ti was multifractal or not. Since each contribution Ti is analog to the transmission of scattering geometry 
of Fig. our present study suggests that each contribution Ti in the incoherent sum of Eq. (22 is not multifractal, 
but monofractal with an exponent Kback — d s = d — 1 in short-range models, i.e. this exponent does not depend on 
the critical multifractal spectrum. 



VI. INFLUENCE OF THE SCATTERING GEOMETRY ON THE TRANSMISSION OFF CRITICALITY 

Up to now, we have only discussed the transmission statistics at criticality. For completeness, we briefly describe 
in this section how the behavior of the transmission depends on the scattering geometry outside criticality. 



A. Influence of the scattering geometry in the localized phase 

In usual short-range models, the localized phase is characterized by exponentially localized wavefunctions, whereas 
in the presence of power-law hoppings, localized wavefunction can only decay with power-law integrable tails. For the 
PRBM, it is moreover expected that the asymptotic decay is actually given exactly by the power-law of Eq. [2] for the 
hopping term defining the model [20| : 



(32) 
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4 5 6 7 8 

In L 

FIG. 7: Statistics of the transmissions T^ 1 ' 1 ', x^ 1 '™' and y£ ls ' m > f or the scattering geometries of Fig. [T] [3] and [5] in the 
localized phase a — 1.4 for b = 0.1 : the measured slopes on this log-log plot are in agreement with Eq. 1331 

As a consequence in the localized phase a > 1, one expects the following typical decays for the transmissions TP' ,1 \ 
T^' m ' and y£ 1B ' m ) f or t ne scattering geometries of Fig. [T] [3] and [5] 

bxT^' 1) (a>l) oc -2alnL 

L — >-oo 

lnT^' m \a > 1) oc -(2a-l)lni 

L — >oo 

ln^ 1B,m) (a > 1) oc -(2a-l)ln£ (33) 

L — >oo 

The decay of T^ 1 ' 1 ^ directly reflects the decay of Eq. [32] of wavefunctions. The presence of many outgoing channels 
change the power by one as a consequence of the integration over the all the sites connected to outgoing channels. 
And the presence of backward channels does not change this exponent. This is in agreement with our numerical 
results shown on Fig. [7] for the case a = 1.4 and b = 0.1. We refer to our previous work [lj| for more details on the 
histogram around the typical value for the one-channel case. 



B. Influence of the scattering geometry in the delocalized phase 

In the delocalized phase, the eigenfunctions are not multifractal, but monofractal with the single value a^e/oc = d 
for the weight I^lC 2 -)!- As a consequence, the typical transmission for the one-channel case of Fig. [I] is expected to 
remain finite as L — > +oo [ItJ (in Eq. [THl the case at yp = d yields Kt yp = 0) 

Tf 1) (a<l) cx TW)>0 (34) 

L — >oo 

The one-channel transmission is thus a good order parameter of the transport properties [l~7j . We refer to our previous 
work [Til for more details on the histogram for this one-channel case. In the presence of many out going wires as in 
Fig. [3l the transmission will thus also remains finite 

T^ m \a <1) oc T( o 1 --)>0 (35) 

L — >oo 

However in the presence of backward scattering channels as in Fig. O one expects that the transmission will now decay 
with L. For instance in the usual short-range Anderson model in dimension d — 3, the many-channel transmission of 
Eq. I5T1 scales as L d ~ 2 , i.e. each contribution in the presence of backward channels scales as 1/L, which represents 
the probability to reach first the out-going surface without returning to the incoming surface for a diffusive particle. 
For the PRBM, our numerical results indicate the typical decay 



lnT| 1B ' m) (a < 1) oc -(2o-l)ln£ (36) 

L — >oo 

which corresponds to the probability to make a direct hopping towards a forward channel without visiting the sites 
connected to the backward channels. 



11 



VII. CONCLUSION 



In this paper, we have studied how the scattering geometry determines the statistics of the Landauer transmission, 
both at criticality and outside criticality. We have presented detailed numerical results for the PRBM model for 
various scattering geometries that interpolate between the two cases that are usually considered, namely the one- 
channel and the many-channel cases. We have found that : (i) in the presence of one isolated incoming wire and 
many outgoing wires, the transmission has the same multifractal statistics as the local density of states of the site 
where the incoming wire arrives; (ii) in the presence of backward scattering channels with respect to the case (i), the 
statistics of the transmission is not multifractal anymore, but becomes monofractal. 



APPENDIX A: REMINDER ON MULTIFRACTALITY OF EIGENFUNCTIONS 



In this Appendix, we recall some useful properties concerning eigenfunction multifractality. The multifractal spec- 
trum f(a) is defined as follows (for more details see for instance the reviews 6, 8]): in a sample of size L d , the number 
Nl{o) of points r where the weight |"0(r)| 2 scales as L~ a behaves as 

Af L {a) oc L /(q) (Al) 

L — >oo 

The typical value at yp corresponds to the maximum value f(a typ ) = d of the multifractal spectrum f(a). The inverse 
participation ratios (I.P.R.s) can be then rewritten as an integral over a 

Y q (L) = [ d d r\^{r)\ 2q ~ / da L f ^ iT« a ~ L~ T ^ (A2) 

J L d J L^oo 

The exponent r(q) can be obtained via a saddle-point calculation in a, and one obtains the Legendre transform 
formula @, Q 

q =/'(<*) 

r(q) =qa- /(a) (A3) 
These scaling behaviors, which concern individual eigenstates ip, can be translated for the local density of states 

PL (E,f)=J2^ E - E nME n (r f )\ 2 (A4) 



as follows : for large L, when the L energy levels become dense, the sum of Eq. IA4I scales as 

p L (E,r)ocL d \i/j E (r)\ 2 (A5) 
and its moments involve the exponents r(q) introduced in Eq. IA2I 

1 



\p L {E,f)]i with A(q)=T(q)-d(q-l) (A6) 

These notions concerning bulk multifractality have been recently extended to surface multifractality [13, |42j 

: the idea is that points near the boundaries are described by another multifractal spectrum fsurfip) which is in 
general not simply related to the bulk multifractal spectrum f(a). 
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